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A density functional theory for colloidal dynamics is presented which includes hydrodynamic 
interactions between the colloidal particles. The theory is applied to the dynamics of colloidal 
particles in an optical trap which switches periodically in time from a stable to unstable confining 
potential. In the absence of hydrodynamic interactions, the resulting density breathing mode, 
exhibits huge oscillations in the trap center which are almost completely damped by hydrodynamic 
interactions. The predicted dynamical density fields are in good agreement with Brownian dynamics 
computer simulations. 
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INTRODUCTION 

The dynamics of mesoscopic colloidal particles dis- 
persed in a molecular solvent is of fundamental impor- 
tance for an understanding of soft matter transport and 
flow properties. The control over the collective colloidal 
dynamics leads to the construction of " smart" materials 
steered by external fields like electro- or magnetorheolog- 
ical fluids [l[ and is essential for applications such as gela- 
tion and aggregation in paints and cosmetics Apart 
from the stochastic Brownian motion of the colloidal par- 
ticles due to their kicks with the solvent molecules, hy- 
drodynamic interactions between colloidal particles aris- 
ing from the induced solvent flow field are getting rele- 
vant for concentrated suspensions. It has been shown by 
experiments [l, [1] , computer simulations [l, Q and the- 
ory 0, H, H, [l^ ^^^^ hydrodynamic interactions can lead 
to qualitatively different behavior in the bulk transport 
properties and in colloidal sedimentation as compared to 
simple Brownian motion valid at very low volume frac- 
tions. 

A full microscopic theory which starts from the col- 
loidal interactions and their hydrodynamic mobility ten- 
sors and predicts the dynamical properties is in principle 
possible by starting from the Smoluchowski picture 
In practice, however, such a predictive theory is ham- 
pered by the many-body nature of the problem and the 
long range of the Oseen mobility tensor which is the lead- 
ing contribution for a colloidal pair. Explicit approaches 
have been worked out in the bulk for short-time and long- 
time diffusion coefficients 0, S B] ^ and for the viscosity 
jl2|. There are also first investigations for colloids near 
walls and on interfaces | ld | l3l| , and for the nonequilib- 
rium structure of colloids [ij] but a general theory for 
an arbitrary and time-dependent inhomogeneous exter- 
nal potential is missing. 

The goal of this letter is twofold: first we construct 
a dynamical density functional theory jl5| which incor- 
porates hydrodynamic interactions. The theory is ex- 



plicitly worked out for hydrodynamic interactions on the 
Rotne-Prager level and generalizes earlier formulations 
[H; [itI . [i8 ] where hydrodynamic interactions were ne- 
glected. The theory makes predictions for an arbitrary 
time-dependent external potential, i.e., for a general in- 
homogeneous nonequilibrium situation. Second, we ap- 
ply the theory to the dynamics of colloidal particles con- 
fined in an unstable optical trap which switches period- 
ically in time from a stable to unstable confining poten- 
tial. This situation can in principle be realized, e.g., by 
combining two laser tweezers or by scannings; around a 
single laser tweezer quickly (l, 0, |2o| . The response to 
this oscillating trap is a time-dependent radial-symmetric 
one-particle density profile which we call - in analogy to 
trapped Bose gases [21] - a driven breathing mode. The 
periodic breathing mode is interesting in itself since it 
may serve as a hydrodynamic transmitter [22| . 

As a result the properties of the breathing mode 
strongly depend on hydrodynamic interactions. For in- 
stance, significant oscillations in the density profile in 
the trap center which built up if no hydrodynamic in- 
teractions are present are completely damped by hydro- 
dynamic interactions. The predictions of the dynami- 
cal density functional theory are in very good agreement 
with Brownian dynamics nonequilibrium computer sim- 
ulations which take hydrodynamic interactions into ac- 
count on the same level as the theory does. The theory 
can in principle be applied to an y in homogeneous sit- 
uation like laser-induced freezing [23]. It may even be 
a possible route to incorporate hydrodynamic interac- 
tions into dynamical approaches like mode-coupling the- 
ory since the latter can be brought into relation with 
dynamical density functional theory [24]. 

The starting point of the derivation of the dynam- 
ical density functional theory including hydrodynamic 
interactions on the two-body level is the Smoluchowski 
equation, i.e., the equation for the time-evolution of the 
full probability density distribution P(r^,t) for N in- 
teracting spherical Brownian particles at positions = 



2 



r 1 , r2 , . . . , r AT and time t (see e.g. ) : 



dt 



N 



(1) 

where k^T is the thermal energy. We assume pair- 
wise additivity for the total potential energy of the 



system, such that U{r^,t) = Y^k^i^e^ti^^k^t) + 
|Z)fcLiE|^/c^2(rfe,rz), where 14xt(r,t) is the one-body 
time-dependent external potential acting on each par- 
ticle and ^2(1*, r') is the pair interaction potential. 
Hydrodynamic interactions are included through the 
configuration-dependent diffusion tensor which we ap- 
proximate on a two particle level: Dij(r^) ^ DolSij + 



'ij Y^^^i - rO + (1 - Sij)u;i2{ri .jj 

Dq denotes the diffusion constant of a single isolated par- 
ticle and 5ij is Kronecker's delta. For a one-component 
suspension of spheres, series expansions of the two ten- 
sors cJi 1 and UJ12 are known, in principle, to arbitrary 
order [25[. By integrating Eq. ([T]) with N J dr2 ... / drAr, 
we obtain the equation for the time-evolution of the one- 
body density piv^t). The resulting equation depends on 
both, the time-dependent two-body and the three-body 
densities. We cast those into a form involving exclu- 
sively the equilibrium Helmholtz free energy functional 
:F[p] = kBT J dr p(r,t)[ln(A3p(r,t)) - 1] + ^exc[p] + 
J dr /)(r, t)T4xt(r, t), with J-exc[p] being the excess con- 
tribution to the free energy functional and A the ther- 
mal de Broglie wavelength, by making use of static DFT 
[26] and the first two members of the Yvon-Born- Green 
(YBG) relations (see, e.g., [l^]). To that end, we iden- 
tify the out-of-equilibrium system at each point in time 
with an equilibrium reference system whose density pro- 
files are identical. The basic assumption now, which also 
underlies the original version of the DDFT [10, [l3 5 is 
to approximate the nonequilibrium two-body and three- 
body densities by those of the reference system with the 
same one-body density. Thus, we obtain our central re- 
sult: 

r-iM£i^)=Vr.L(M)v ^^[^] . 



Here, 



dt 



5p{r,t) 
5:F[p] 



|drV(2)(r,r',t)a;i2(r-r')-V: ^^^"^ 



5p{v',t) ' 



(2) 



with the mobility constant F for which the Einstein re- 
lation gives Dq/T = (/cbT)~^. Eq. ([2]) has the form 
of a continuity equation with the current density j = 
ji+j2+j3 given by the terms in the curly brackets. The 
current density ji is proportional to the thermodynamic 



6pir,t) 



and persists when hydrodynamic inter- 



force V, 

actions are neglected pj|, Il8|. j2 and js are additional 



current densities which occur due to the solvent medi- 
ated hydrodynamic interactions. J2 describes the current 
density stemming from the refiection of the solvent fiow 
induced by the thermodynamic force at position r on the 
surrounding particles, js, on the other hand, is the cur- 
rent density at position r due to the solvent fiow induced 
by the thermodynamic force at position r^ 

Finally, we close the above relation Eq. ([2]), which 
still depends on the nonequilibrium two-body density 
p^'^\r^r\t). Within our approximation - the two-body 
density is assumed to be identical to the equilibrium one 
of the reference system - it is given at every point in time 
by the exact generalized Ornstein-Zernike equation [2^: 

p(2)(r,r',t) = p{r,t)p{r',t) (l + (kBT)-^ j£^i^) 4 

p{r', t) J dr" I (r, r", t) - p{r, t)p{r", t)) 



{kBT) 



-1 S^J^e^cM 

Sp{r',t)Sp{r 



(3) 



This implicit equation for the two-body density of 
the inhomogeneous system ma y b e reasonably approx- 
imated by its bulk value [H, [2^, i.e.: p^'^\r,r',t) ^ 
p{r^t)p{r' ^t)g{\r — r'|,p), where ^(|r — r'|,p) is the pair 
correlation function for a homogeneous system at an ap- 
propriately averaged density p. For a hard sphere fiuid, 
an analytic expression for the pair correlation function is 
available based on the Percus- Yevick equation [s^ . 

We use the method presented here to investigate the 
time-evolution of the one-body density of a confined clus- 
ter of = 100 monodisperse hard spherical particles of 
diameter a, which serves as the unit of length hence- 
forth. An appropriate time scale is tb = ct^/I^o, and the 
energy unit is /cbT. The particles are trapped in a soft 
spherical cavity which switches from a stable to an un- 
stable shape periodically in time. The confining external 
potential only acts on the colloidal particles. Therefore 
the solvent is treated as an unbounded fiuid. The total 
external potential is modeled as 



F2Cos(27rt/T)( — ) , (4) 



where r = |r|, i?i = 4<j and Vi — 10/cbT are the length 
scale and the strength of an outer fixed cavity and R2 = cf 
and V2 = k^T are the length scale and strength of an 
inner part, which oscillates in time with a period r = 
O.Stb- a sketch of the setup is shown in Fig.[Tl Due to the 
spherical symmetry, the density profile p(r, t) depends 
only on the radial position coordinate r. 

For the hard sphere excess density functional ^exc[p] 
Rosenfeld's fundamental measure theory |32| was used 
which provides a very reliable approximation in equi- 
librium [3^. The distinct hydrodynamic tensor cJi ofr) 
is approximated by the Rotne-Prager expression [3J] 
^ou;i2(r) = + rr] + - 3ff] + 0[{^y] 




FIG. 1: (Color online) Sketch of the confined system. The 
external potential models an optical trap Vext{r,t) which 
changes its central shape from stable to unstable within a 
time period r. The trap confines N colloidal hard spheres of 
diameter a shown as black circles. Additionally, typical 3d 
simulation snapshots are shown. The left hand side shows an 
initial stable configuration for t — and the right hand side 
shows an unstable situation at t = 2.75tb for case (N) [3l| . 



while the self term u;ii(r) whose leading order term is 
0((cr/r)^) is neglected. Here, r = r/|r| denotes the unit 
vector, rr is the dyadic product, and 1 is the unit ma- 
trix. Thus, on this level of approximation, we incor- 
porate all solvent mediated interactions up to order of 
0{{a/r)^). The pair correlation g{\T — t'\^ p) is calcu- 
lated at each time step at the average density of the 
system p{t) = 1/Rm.^{t) J^'^^''^'^ drp{r,t), where i^^ax 
is defined by Ve^t{r = i^max(^)) = IO/crT. 

The results are tested against Brownian dynamics sim- 
ulations [sHI performed on the same level of accuracy of 
the diffusion tensor, in which the hard interaction is ap- 
proximated by a slightly softened one: 



V2{r) 
kBT 



(7) 



X 48 







if r < 
else 



(5) 



In all simulations we chose a finite simulation time step of 
At = IO^^tb- In order to obtain the time-dependent den- 
sity /)(r, t) we perform a large number of A^run = 10^ inde- 
pendent runs with different initial configurations sampled 
from a situation with a static external potential, i.e., Eq. 
(m at t = 0. Additionally, the densities are compared to 
those obtained by standard DDFT where hydrodynamic 
interactions are ignored, i.e., = (jOi2 = 0. Henceforth, 
we label the situation including hydrodynamic interac- 
tions (H) and the situation where they are neglected (N), 
respectively. The initial density profile is the equilibrium 
density profile for 14xt(^, ^ = 0) calculated from static 
density functional theory. Typical simulation snapshots 
are shown in Fig. [H 

The resulting steady-state of the dynamical density 
profiles, i.e. the driven breathing mode after initial re- 
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FIG. 2: (Color online) Steady-state DDFT (solid curves) and 
BD (noisy curves) results for the time dependent density pro- 
file p{r,t). In Fig. (a) and (b) hydrodynamic interactions are 
taken into account while in (c) and (d) they are neglected, (a) 
and (c) correspond to the expanding half period and (b) and 
(d) to the compressing half period, respectively [31]. The pro- 
files correspond to the following time sequence: to = 2.5tb, 
ti = 2.6rB, t2 = 2.7rB, ts = 2.75rB in (a) and (c), and 
t4 = 2.75TB, ts = 2.85rB, ts = 2.9tb, and te = S.Otb. 



laxation, is depicted in Fig. [2l First of all, theory and 
simulation results are in very good agreement for both 
situations (H) and (N) but we observe distinct qualita- 
tive differences in the breathing mode: Hydrodynamic 
interactions tend to damp the density response consid- 
erably. For neglected hydrodynamic interactions there 
are huge density oscillations, in particular at the trap 
center but also at the trap boundaries which are signifi- 
cantly smaller for hydrodynamic interactions. This result 
is not obvious as hydrodynamic interactions tend to ac- 
celerate neighboring particle which are driven into the 
same direction. The damping effects seen here is caused 
by the overall motion of the breathing mode which hin- 
ders collective streaming due to the counter motion in 
the opposed part of the trap. 

In order to analyze the relaxational behavior towards 
the steady state, we introduce the second moment of the 
breathing mode m2{t) = J dr r'^p{r^t). It is shown in 
Fig. [3] in both cases (H) and (N) for DDFT and Brow- 
nian Dynamics. Clearly, the dynamic evolution of the 
second moment is strongly damped by hydrodynamic 
interactions as revealed by the much slower oscillation 
amplitude. On top of that the relaxation time towards 
the steady state breathing mode is considerably larger 
for hydrodynamic interactions as compared to the sim- 
ple Brownian case where the relaxation is almost instan- 
taneous. The second moment of the breathing mode is 
slightly off-phase with respect to the driving external po- 
tential ~ cos(27rt/r), and the hydrodynamic interactions 
lead to a stronger phase-shifting. 



4 




I \ \ \ \ \ \ \ \ \ \ \ I 

0.5 1 1.5 2 2.5 3 



FIG. 3: (Color online) Second moment of the breathing mode, 
7712 (t), versus time t. DDFT (solid curves) and BD (dashed 
curves) results for hydrodynamic interactions taken into ac- 
count (H) and being neglected (N). 



In conclusion, we have proposed a dynamical density 
functional theory which includes hydrodynamic interac- 
tions between the colloidal particles and applied it to 
access the driven breathing mode in oscillating optical 
traps. The theory was confirmed by Brownian dynamics 
computer simulations. Hydrodynamic interactions were 
found to damp the response to the driving trap, to in- 
crease the relaxation time towards the steady state and to 
increase the phase shift. These predictions can in princi- 
ple be tested by real-space experiments on confined col- 
loidal particles. For charged suspensions, the strength 
of the hydrodynamic interactions can be systematically 
tuned in the experiments by varying the colloid charge 
which governs size of the interactions relative to the phys- 
ical core which governs the hydrodynamic interactions. 
It would be interesting to generalize the theory further 
to orientational degrees of freedom and to self-propelled 
colloidal particles representing microswimmers where hy- 
drodynamic interactions play a key role. 
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